in 



Faster exponential-time algorithms in graphs of bounded average 

degree* 

Marek Cygant Marcin Pilipczuk-'- 

February 18, 2013 

1. 

. We first show that the Traveling Salesman Problem in an n-vertex graph with average degree bounded 

["T | . by d can be solved in and exponential space for a constant Ed depending only on d. 

Thus, we generalize the recent results of Bjorklund et al. [TALG 2012] on graphs of bounded degree. 

Then, we move to the problem of counting perfect matchings in a graph. We first present a simple 
algorithm for counting perfect matchings in an n-vertex graph in 0*(2 n//2 ) time and polynomial space; 
our algorithm matches the complexity bounds of the algorithm of Bjorklund [SODA 2012], but relies 
^0 ■ on inclusion-exclusion principle instead of algebraic transformations. Building upon this result, we show 

that the number of perfect matchings in an n-vertex graph with average degree bounded by d can be 
computed in O* (2 (1 ^ £2d ^"^ 2 ) time and exponential space, where E2d is the constant obtained by us for 
. the Traveling Salesman Problem in graphs of average degree at most 2d. 

Moreover we obtain a simple algorithm that counts the number of perfect matchings in an n-vertex 
bipartite graph of average degree at most d in £)*(2' 1 ~ 1/ '' 3 ' 55d - ) - )n ' /2 ) time, improving and simplifying the 
recent result of Izumi and Wadayama [FOCS 2012]. 

CO 

^ ■ 1 Introduction 

CO 

Improving upon the 50-years old 0*(2™)-time dynamic programming algorithms for the Traveling Salesman 
Problem by Bellman [I] and Held and Karp [7. is a major open problem in the field of exponential-time 
, algorithms [14] . A similar situation appears when we want to count perfect matchings in the graph: a half- 

century old 0*(2"/ 2 )-time algorithm of Ryser for bipartite graphs [H] has only recently been transferred to 
arbitrary graphs [3J, and breaking these time complexity barriers seems like a very challenging task. 

From a broader perspective, improving upon a trivial brute-force or a simple dynamic programming 
algorithm is one of the main goals the field of exponential-time algorithms. Although the last few years 
brought a number of positive results in that direction, most notably the 0*(1.66") randomized algorithm for 
finding a Hamiltonian cycle in an undirected graph it is conjectured (the so-called Strong Exponential 
Time Hypothesis [5]) that the problem of satisfying a general CNF-SAT formulae does not admit any 
exponentially better algorithm than the trivial brute-force one. A number of lower bounds were proven 
using this assumption [6j [101 E] ■ 

In 2008 Bjorklund et al. |S] observed that the classical dynamic programming algorithm for TSP can be 
trimmed to running time 0*(2( 1_eA ) n ) in graphs of maximum degree A. The cost of this improvement is the 
use of exponential space, as we can no longer easily translate the dynamic programming algorithm into an 
inclusion-exclusion formula. The ideas from [5] were also applied to the Fast Subset Convolution algorithm, 
yielding a similar improvements for the problem of computing the chromatic number in graphs of bounded 
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degree [3]. In this work, we investigate the class of graphs of bounded average degree, a significantly broader 
graph class than this of bounded maximum degree. 

In the first part of our paper we generalize the results of [5], 

Theorem 1.1. For every d > 1 there exists a constant Sd > such that, given an n-vertex graph G of 
average degree bounded by d, in O* (2( 1 ~ e< *) n ) time and exponential space one can find in G a smallest weight 
Hamiltonian cycle. 

We note that in Theorem 1 1 . 1 1 the constant e c i depends on d in doubly-exponential manner, which is worse 
than the single-exponential behaviour of [5 in graphs of bounded degree. 

The proof of Theorem 11.11 follows the same general approach as the results of [5] — we want to limit 
the number of states of the classical dynamic programming algorithm for TSP — but, in order to deal with 
graphs of bounded average degree, we need to introduce new concepts and tools. Recall that, by standard 
averaging argument, if the average degree of an n-vertex graph G is bounded by d, for any D > d there are 
at most dn/D vertices of degree at least D. However, it turns out that this bound cannot be tight for a large 
number of values of D at the same time. This simple observation lies at the heart of the proof of Theorem 
II. li as we may afford more expensive branching on vertices of degree more than D provided that there are 
significantly less than dn/D of them. 

In the second part, we move to the problem of counting perfect matchings in an n-vertex graph. Wc 
start with an observation that this problem can be reduced to a problem of counting some special types 
of cycle covers, which, in turn, can be done in C*(2™ //2 )-time and polynomial space using the inclusion- 
exclusion principle (see Section 15. ip . Note that an algorithm matching this bound in general graphs has 
been discovered only last year [3], in contrast to the 50-years old algorithm of Ryser for bipartite graphs. 
Thus, we obtain a new proof of the main result of [3j, using the inclusion-exclusion principle instead of 
advanced algebraic transformations. 

Once we develop our inclusion-exclusion-based algorithm for counting perfect matchings, we may turn it 
into a dynamic programming algorithm and apply the ideas of Theorem 11.11 obtaining the following. 

Theorem 1.2. Given an n-vertex graph G of average degree bounded by d, in O* ' (2^ 1 ~ e2d ' n ^ 2 ) time and 
exponential space one can count the number of perfect matchings in G where S2d is the constant given by 
Theorem \l.l\ for graphs of average degree at most 2d. 

To the best of our knowledge, this is the first result that breaks the 2™/ 2 -barrier for counting perfect 
matchings in not necessarily bipartite graphs of bounded (average) degree. 

When bipartite graphs are concerned, the classical algorithm of Ryser [T2] has been improved for graphs 
of bounded average degree first by Servedio and Wan [13] and, very recently, by Izumi and Wadayama [5]. 
Our last result is the following theorem. 

Theorem 1.3. Given an n-vertex bipartite graph G of average degree bounded by d, in 0*(2^ 1 ^ 3 - 55d ^ n / 2 ) 
time and exponential space one can count the number of perfect matchings in G. 

Hence, we improve the running time of [51 113) in terms of the dependency on d. We would like to 
emphasise that our proof of Theorem 1 1.3 1 is elementary and docs not need the advanced techniques of coding 
theory used in [9]. 

Organization of the paper Section [2] contains preliminaries. Next, in Section [3] we prove the main 
technical tool, that is Lemma [3^11 used in the proofs of Theorem 11.11 and Theorem 1 1.21 In Section 2] we prove 
Theorem 11.11 while in Section 15.11 we first show an inclusion-exclusion based algorithm for counting perfect 
matchings, which is later modified in Section 15.21 to fit the bounded average degree framework and prove 
Theorem 11.21 Finally, Section [5] contains a simple dynamic programming algorithm, proving Theorem II. 31 

We would like to note that both Section 15.11 and Section [5] are self-contained and do not rely on other 
sections (in particular do not depend on Lemma |3.4[) . 
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2 Preliminaries 



We use standard (multi)graph notation. For a graph G = (V, E) and a vertex v G V the neighbourhood of 
v is defined as Ng(v) — {u : uv G £?} \ {w} and the closed neighbourhood of v as A^v] = Ng{v) U {w}. The 
degree of v € V is denoted deg G (w) and equals the number of end-points of edges incident to v. In particular 
a self-loop contributes 2 to the degree of a vertex. We omit the subscript if the graph G is clear from the 
context. The average degree of an n-vertex graph G = (V,E) is defined as — YlveV deg(v) = 2\E\/n. A 
cycle cover in a multigraph G = (V, E) is a subset of edges C C E, where each vertex is of degree exactly 
two if G is undirected or each vertex has exactly one outgoing and one ingoing arc, if G is directed. Note 
that this definition allows a cycle cover to contain cycles of length 1, i.e. self-loops, as well as taking two 
different parallel edges as length 2 cycle (but does not allow using twice the same edge). 

For a graph G — (V, E) by Vd og =c, Vd C g> c , Vd og > c let us denote the subsets of vertices of degree equal to 
c, greater than c and at least c respectively. 

We also need the following well-known bounds. 

Lemma 2.1. For any n,k > 1 it holds that 

(;)*(?)*■ 

Lemma 2.2. For any n > 1, it holds that H n —i > lnn 7 where H n — X)"=i J- 

Proof. It is well-known that linijj^oo H n — Inn = 7 where 7 > 0.577 is the Euler-Mascheroni constant and 
the sequence H n — \nn is decreasing. Therefore H n -\ = H n — — > Inn + 7 — i, hence the lemma is proven 
for n > 2 as 7 > h. For n=l, note that H n _i = Inn = 0. □ 



3 Properties of bounded average degree graphs 

This section contains technical results concerning bounded average degree graphs. In particular we prove 
Lemma 13.41 which is needed to get the claimed running times in Theorems 11.11 and 11.21 However, as the 
proofs of this section arc not needed to understand the algorithms in further sections the reader may decide 
to see only Definition 13.31 and the statement of Lemma [ 



Lemma 3.1. Given an n-vertex graph G = (V, E) of average degree at most d and maximum degree at most 
D one can in polynomial time find a set A containing |~ 2+ " rfg ] vertices of degree at most 2d, where for each 
x,y e A, x ^ y we have Ng[x\ H Nc[y] = 0- 

Proof. Note that |Vd eg <2d| > n/2. We apply the following procedure. Initially we set A := and all the 
vertices are unmarked. Next, as long as there exists an unmarked vertex x in Vd eg <2d, we add x to A and mark 
all the vertices 7Vg[7Vg[2;]]. Since the set Ng[Ng[x]} contains at most 1 + 2d+ 2d(D — 1) = 1 + 2dD vertices, 
at the end of the process we have |A| > 2 +ldD • Clearly this routine can be implemented in polynomial 
time. □ 

Lemma 3.2. For any a > and an n-vertex graph G = (V, E) of average degree at most d there exists 
D < e a such that \V dcf r >D \ < 

Proof. By standard counting arguments we have 

00 00 

^2 \Vde S >i\ = ^2i\V dcg=l \ < nd. 

i=a i=o 

For the sake of contradiction assume that |Vdeg>i| > , for each i < e a . Then 

00 L e °J L e °J n ^ 

^ |Vdeg>t| > ^ l^dcg>i| > V" l/i = i?Le Q J ^ nd > 

i=0 i=l i=l 
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where the last inequality follows from Lemma \2. 2 1 □ 

In the following definition we capture the superset of the sets used in the dynamic programming algorithms 
of Theorems O and IL"2l 



Definition 3.3. For an undirected graph G = (V, E) and two vertices s,t £ V by deg2sets(G, s, t) we define 
the set of all subsets X C\ V \ {s, t}, for which there exists a set of edges F C E such that: 

• deg F {v) = for each v G V \ (X U {s,t}) : 

• deg F (v) = 2 for each v £ X, 

• &eg F (v) < 1 for v G {s,t}. 

Lemma 3.4. For every d > 1 there exists a constant Ed > 0, such that for an n-vertex graph G = (V, E) of 
average degree at most d for any s,t <E V the cardinality of deg2sets(G, s, t) is at most O* (2^~ ed ^ n ) . 

Proof. Use Lcmma l3~l2"l with a — e cd for some sufficiently large universal constant c (it suffices to take c = 20). 
Hence we can find an integer D < e a = e e such that there are at most vertices of degree greater than 
D in G. 

Let D' = max(2d,D) and H = G[Vd og <D']- Moreover let Y = V dcg>D > and recall \Y\ < as D' > D 
and Y C Vdeg>£>- Note that H contains at least n/2 vertices and has average degree upper bounded by d. 
By Lemma |3~T1 there exists a set A C V(H) of \n/(4 + 8dD')~\ vertices having disjoint closed neighbourhoods 
in H . Note that, since d > 1 and D' > 2d: 



\A\ 



+ 8dD 

If n < 8 ?_ , n — 0(1) and the claim is trivial. Otherwise: 



I li IV IV / -t \ 

~ 4 + 8dD' ~ 2dD' + MD' ~ IQdD' ' ^ ^ 



< 8^ + 1< 2^ (2) 



4 + 8dD 

Moreover, as d > 1 and D' = max(2<i, D) < 2eLD, for sufficiently large c we have: 



1 1 ~ aD ~ 2MD' e cd 2MD 1 ~ 2 ' w 



Consider an arbitrary set X G deg2sets(G, s, t), and a corresponding set F C £ from Definition 13.31 
Define Zx as the set of vertices x G X n V(£f) such that Njj(x) D X = 0. Note that F is a set of paths and 
cycles, where each vertex of Zx is of degree two, hence F contains at least 2\Zx\ edges between Zx and Y, 
as any path/cycle of F visiting a vertex of Zx has to enter from Y and leave to Y. Hence by the upper 
bound of 2 on the degrees in F we have \Zx\ < \Y\. 

For each x G A \ (Z x U {M}) we have that AT^ar] n X ^ {x} and |iV H [x]| < 2rf + 1. By definition, if 
x E AD Zx, we have Njj[x] D X = {x}. Therefore, for fixed v and A n there are at most 

j2d+l _ I \ l^\(ZxU{ S ,t})| / v |AOZx| /,2<H-1 ,\W 

2" | - „ , , ) l-s^rl <2"+ 2 ' ' 



22d+i / I 2 2d+1 7 — V 2 2d+1 

choices for X G deg2sets(G, s). 

Moreover, there are at most X)|=o ('f ) — n (jyJ) choices for Zx D A. Thus 



| ( I, A 2s, ( s(G.>,/) <2^.( " 2M+1 J -r^j) (4) 
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Let us now estimate ( £ ) by Lemma |2~T1 Since L> < D', \Y\ < and by © and ©: 
f\A\\ fe\A\\ lYl ( n aD\^ / a \£* 



By the standard inequality I — x < e x we have that 

(2 2d+l _ 1)/2 2d+l = (1 _ ^22^+1) < e -l/2*H-* _ (6) 

Using ([J), (|5|) and ([6]) we obtain that 

(\A\\ f 2 2d+1 - 1 \ |A|/2 /ndlna n 

^ 2 2d +! J " 6XP ^^D - ~ 20dD'2 2d + 1 

Plugging in a = e cd and using the fact that e 10d > 40d 2 for d > 1 we obtain: 

/2 2d+1 - i\ |A|/2 / ncd n 



\Y\J V 2 2d +! y - 6XP Ve( c - 10 ) d 20d-2di^ 2(W2 2d + 1 
Since D' = max(2d,L>) < 2dL>' and e 4d > 2 2d+1 as d > 1, we get 

/L4|\ /2^+i^l\ |A|/2 / n / c _ u 

Vl^ly V 22d+1 / " eXP V20d^'2 2 <i+i U( c - 14 ) d //' 
Finally, for sufficiently large c, as d > 1, we have c < e( c ~ 14 ' d and 

L4|\ /2 2d+1 - 1 



,2(1+1 1 \ l^l/ 2 



Consequently plugging ([7]) into (@| and using (JTJ) and © we obtain: 

- 2 2d+l _ l^l/ 2 



< 1. (7) 



|deg2sets(G, s,t)\ < n2 



n+2 



2 2(i+l 

< n2 n+2 exp (—5^ ) 

f\ 2 2d+i . 20dD'/ 

< n2"+ 2 exp / 



2 2d +! • 20d ■ e eCd - 

This concludes the proof of the lemma. Note that the dependency on d in the final constant Ed is doubly- 
exponential. □ 



4 Algorithm for TSP 

To prove Theorem it suffices to solve in O* (2^~ ed ^ n ) time the following problem. We are given an 
undirected n-vertex graph G = (V, E) of average degree at most d, vertices a, b € V and a weight function 
c : E — > M+. We are to find the cheapest Hamiltonian path between a and b in G, or verify that no 
Hamiltonian a6-path exists. 

We solve the problem by the standard dynamic programing approach. That is for each a G X C V and 
v E X we compute t[X][t;], which is the cost of the cheapest path from a to v with the vertex set X. The 
entry i[V][6] is the answer to our problem. Note that it is enough to consider only such pairs (X, v), for 
which there exists an aw-path with the vertex set X . 

We first set i[{a}][a] = 0. Then iteratively, for each i = 1, 2, . . . , n — 1, for each u G V, for each X C V 
such that \X\ = i, a,u € X and t[JC][tt] is defined, for each edge uv G E where v $ X, if i[X U is 
undefined or tLY U > t[X][u] + c(m>), we set t[X U {«}][«] = t[X][u] + c(uv). 

Finally, note that if t[JC][v] is defined then X \ {a, v} € deg2sets(G, a, u). Hence, the complexity of 
the above algorithm is within a polynomial factor from X^ev |deg2sets(G, a, v)\, which is bounded by 
0*(2 (1 - e ^ n ) by Lemma [3H 
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5 Counting Perfect Matchings 



In this section we design algorithms counting the number of perfect matchings in a given graph. First, in 
Section 15.11 we show an inclusion-exclusion based algorithm, which given an n- vertex graph computes the 
number of its perfect matchings in 0*(2"/ 2 ) time and polynomial space. This matches the time and space 
bounds of the algorithm of Bjorkhmd [3,. Next, in Section I5TH1 we show how the algorithm from Section I5~T1 
can be reformulated as a dynamic programming routine (using exponential space), which together with 
Lemma 13.41 will imply the running time claimed in Theorem 11.21 



5.1 Inclusion-exclusion based algorithm 

In the following theorem we show an algorithm computing the number of perfect matchings of an undirected 
graph in 0*(2™/ 2 ) time and polynomial space, thus matching the time and space complexity of the algorithm 
by Bjorklund [3]. 

Theorem 5.1. Given an n-vertex graph G — (V,E) in 0*{2 n / 2 ) time an polynomial space one can count 
the number of perfect matchings in G. 

Proof. Clearly we can assume that n is even. Consider the edges of G being black and let V = {vo, . . . , v n ~i}. 
Now we add to the graph a perfect matching of red edges Er = {v2iV2i+i ■ < i < n/2} obtaining 
a multigraph G' . Observe that for any perfect matching MCE the multiset M U En is a cycle cover 
(potentially with 2-cycles), where all the cycles are alternating - that is when we traverse each cycle of 
M U En, the colors alternate (in particular, they have even length). Moreover, for any cycle cover Y of G" 
composed of alternating cycles the set Y \ En is a perfect matching in G. This leads us to the following 
observation. 

Observation 5.2. The number of perfect matchings in G equals the number of cycle covers in G' where 
each cycle is alternating. 

Now we create a directed multigraph graph G" with arcs labeled with elements of L = {£q, . . . ,£ n /2-i}i 
having n vertices and 2m arcs, where m = \E\ is the number of black edges of G". Let {v'q , . . . , be the 

set of vertices of the graph G" . For each black edge u a u& of G" we add to G" two following arcs: 

• «biX') labeled £ la/2 ], 

• and «' ffil , v") labeled £y b/2 \- 

By © we denote the XOR operation, that is, for any < x < n the vertex v x( ^\ is the other endpoint of the 
red edge of G 1 incident to v x . 

Observation 5.3. The number of cycle covers in G' where each cycle is alternating equals the number of 
sets of cycles in G" of total length n/2, where each label li (for < i < n/2) is used exactly once. 

We are going to compute the of sets of cycles in G" where each label is used exactly once using the 
inclusion-exclusion principle. 

For a vertex i>" of G", we say that a closed walk C is v"-nice if C visits v" exactly once and does not 
visit any vertex v'l for b < a. A closed walk is nice if it is w"-nice for some w"; note that, in this case, the 
vertex v 1 ^ is defined uniquely. For a positive integer r let us define the universe Q r as the set of r-tuples, 
where each of the r coordinates contains a nice closed walk in G" and the total length of all the walks 
equals n/2. For < i < n/2 let A T: i C fl r be the set of r-tuples, where at least one walk contains an 
arc labeled £j. Note that by the observations we made so far the number of perfect matchings in G equals 
Si<r<n/2 I no<i<n/2 1 1 r ^ i as the tuples in Q r are ordered and in any tuple of Hckioi/2 -^r.i a ^ wa lks are 
pairwisc different. Therefore from now on we assume r to be fixed. By the inclusion-exclusion principle 



n • 

0<i<n/2 



7C{0,...,n/2-l} 



G 



hence to prove the theorem it is enough to compute the value | P| ieJ (Sl T . \ A^j) | for a given / C {0, . . . , n/2 — 1} 
in polynomial time. Let G'[ be the graph G" with all the arcs with a label from Lj — \li : i G /} 
removed. Let p a j be the number of v'^-nice closed walks in G" of length j. Note that the value p a j can 
be computed in polynomial time by standard dynamic programming algorithm, filling in a table ip[fr][i], 
a < b < n/2, < i < j, where i P [&][j] is the number of walks W from v" to v' b ' in G" of length i that visit 
only once and does not visit any vertex i>" for c < a. 

Finally, having the values p a j is enough to compute | P) ie 7(^-r \ -A*) I by the standard knapsack type 
dynamic programming. That is, we fill in a table < q < r, <i < n/2, where t [g] [i] is the number 

of g-tuplcs of nice closed walks in G'[ of total length i. 

□ 

5.2 Dynamic programming based algorithm 

To prove Theorem 1 1.2 1 we want to reformulate the algorithm from Section [5.11 to use dynamic programming 
instead of the inclusion exclusion principle. This causes the space complexity to be exponential, however it 
will allow us to use Lemma 13.41 to obtain an improved running time for bounded average degree graphs. 

Assume that we are given an n- vertex undirected graph G ~ (V, E), where n is even, and we are to count 
the number of perfect matchings in G. We are going to construct an undirected multigraph G' having only 
n/2 vertices, where the edges of G' will be labeled with unordered pairs of vertices of G', i.e. with edges of 

G. As the set of vertices of G' = (V, E 1 ) we take V — {v' , . . . , v' n ^_ 1 }. For each edge v a Vb of G we add to 
G' exactly one edge: v \ a /2\ v [b/2] labeled with {v a , Vb}- For an edge e' G E' by £(e') let us denote the label of 
e'. Note that G' may contain self- loops and parallel edges. Observe that if the graph G is of average degree 
d, then the graph G' is of average degree 2d. 

In what follows we count the number of particular cycle covers of G', where we use the labels of edges 
to make sure that a cycle going through a vertex v[ G V never uses two edges of G' corresponding to two 
edges of G incident to the same vertex. 

Lemma 5.4. The number of perfect matchings in G equals the number of cycle covers C C E' of G' , where 

Proof. We show a bijection between perfect matchings in G and cycle covers G of G' satisfying the condition 
{J eeC t{e) = V. 

Let M be a perfect matching in G. As f(M) we define f{M) — { v 'y a / 2 \ v '\b/2\ '■ VaVb e ^O- Note that 
f(M) is a cycle cover and moreover U e e/(M) ^( e ) = V ■ ^ n the reverse direction, for a cycle cover G C E' 
of G', consider a set of edges h(C) defined as h(C) — {i{e) : e G G}. Clearly the condition \J eeC £(e) = V 
implies that h(C) is a perfect matching, and moreover h = □ 

Observe, that if a cycle cover G C E' of G' does not satisfy Ueec ^( e ) = ^> then there is a vertex v[ G V, 
such that the two edges of G incident to v\ do not have disjoint labels. Intuitively this means we are able to 
verify the condition U e ec^( e ) = ^ locally, which is enough to derive the following dynamic programming 
routine. 

Lemma 5.5. Once can compute the number of cycle covers C of G' satisfying [J eeC £{e) — V in 
0*(J2 S t£V |deg2sets(G', s, t)\) time and space. 

Proof. An ordered r-cycle cover of a graph H is a tuple of r cycles in H , whose union is a cycle cover of 

H . As each cycle cover of H that contains exactly r cycles can be ordered into exactly r! different ordered 
r-cycle covers, it is sufficient to count, for any 1 < r < n/2, the number of ordered r-cycle covers G in G' 
such that each two edges in G have disjoint labels. In the rest of the proof, we focus on one fixed value of r. 

For < q < r and X C V' as t[q][X] let us define the number of ordered g-cycle covers in G'[X] where 
each two edges have disjoint labels; note that t[r][V'] is exactly the value we need. Moreover for < q < r, 
X C V , v' a ,v' b G X, a < b and x G {w2b,i'2h+i} as t 2 [<z][A'][wQ][v[,][a;] we define the number of pairs (G, P) 
where 
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• C is a ordered g-cycle cover of G'[Y] for some Y C X \ {v' a , v' b }; 

• P is a i/ i^-path with vertex set X \ y that does not contain any vertex v' c with c < a; 

• any two edges of C U P have disjoint labels; 

• the label of the edge of P incident to v' a contains v 2a \ 

• the label of the edge of P incident to v' b contains x. 

Note that we have the following border values: t[0] [0] = 1 and t[0][X] = for X ^ 0. 

Consider an entry t 2 [q] [X] [v' a ] [v' b ] [x] , and let (C, P) be one of the pairs counted in it. We have two cases: 
either P is of length 1 or longer. The number of pairs (C, P) in the first case equals t[q] [X \ {v' a , v' b }] ■ \{v' a v' b G 
E' : £(v' a v b ) = {v2a, x}}\- In the second case, let v' c v' b be the last edge of P; note that c > a by the assumptions 
on P. The label of v' c v' b equals {v2 C ,x} or {v 2c +i,x}. Thus, the number of elements (C,P) in the second 
case equals E„/ex\K,«J} E y e{v 2c ,v 2c+1 } Ho\{ x \ KMKKlfo © 1] • \WA e E' : £(v' c v' b ) = {y,x}}\, where 
for y — v' r we define y © 1 = u^i- 

Let us now move to the entry and let C be an ordered g-cycle cover in G'[X]. Again, there are 

two cases: either the last cycle of C (henceforth denoted W) is of length 1 or longer. The number of the 
elements C of the first type equals J2 V ' ex % — \ K}] • |{ w l w a e ^ n tne second case, let v' a be 

the lowest-numbered vertex on W and let e = v' a v' b be the edge of W where v 2a+ i € £(e). Note that both v' a 
and e are defined uniquely; moreover, a < b and no vertex u£ with c < a belongs to W. Thus the number 
of elements C of the second type equals ^2 v / atV > b ^x,a<b^2xe{v 2 b,v2b+i} 

t 2 [q-l][X}[v' a }[v' b }[x®l}-\{v' a v' b eE' : 

Z( v 'a v 'b) = {V2a+l,x}}\. 

So far we have given recursive formulas, that allow computing the entries of the tables t and t 2 . However 
the values t% [q] [X] [v' a ] [v' b ] [x] for X {J s t v , deg2sets(G", s, t) arc equal to zero. The last step of 

the proof is to show how to perform the dynamic programming computation in a time complexity within 
a polynomial factor from the number of non-zero entries of the table. We do that in a bottom-up manner, 
that is iteratively, for each q = 1, 2, . . . , r, for each i = 1,2, . . . , n, we want to compute the values of non-zero 
entries f[g][X] for all sets X of cardinality i and then compute the values of non-zero entries t 2 [q] [X] [*][*] [*] 
for all sets X of cardinality i. Having the non-zero entries for the pairs (q',i') where q' < q, i < i one 
can compute the list of non-zero entries for \X\ = i by investigating to which recursive formulas the 

non-zero entries for (q',i') contribute to. Analogously having the non-zero entries for the pairs {q 1 ' where 
q 1 < q, i' < i we generate the non-zero entries t<i [q] [X] [*][*] [*] for \X\ = i, which finishes the proof of the 
lemma. 

□ 

Theorem 11.21 follows directly from the Lemma 13.41 together with Lemma 15.51 

6 Counting Perfect Matchings in Bipartite Graphs 

In this section we prove Theorem 11.31 i.e. show an algorithm counting the number of perfect matchings in 
bipartite graphs of average degree d in 

0*(2(i-l/(3.55d)W2) time; 

improving and simplifying the algorithm 

of Izumi and Wadayama [9] . 

Let G = (V = AUB, E) be a bipartite graph, where \A\ = \B\ = n/2, and denote k = n/2. Note that we 
may assume that each vertex in G is of degree at least 2, as an isolated vertex causes no perfect matching to 
exist, while a vertex of degree 1 has to be matched to its only neighbour, hence we can reduce our instance 
in that case. Therefore we assume d > 2. 

Let Bo C B be a subset containing [k/(ad)\ vertices of smallest degree in B, where a > 2 is a constant 
to be determined later. Moreover let A = N(B ) and observe that \A \ < k/a, as vertices of B are of 
average degree at most d. We order vertices of A, i.e. denote A — {a 1; . . . , a^}, so that vertices of A \ A 
appear before vertices of A . In particular for any 1 < i < k(l — l/a) we have N(a.i) n B = 0. 
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Consider the following standard dynamic programming approach. For X C B define t[X] as the number 
of perfect matchings in the subgraph of G induced by {ai, . . . , am} UX. Having this definition the number 
of perfect matchings in G equals T[B]. Observe that the following recursive formula allows to compute the 
entries of the table t, where we sum over the vertex matched to a,\x\'- 

t[x\= ]T t[x\{v}\, 

where i[0] is defined as 1. 

Let us upper bound the number of sets X, for which t[X] is non-zero. If \X\ < (1 — l/a)k and t[X] > 0, 
then X n Bq = 0, as otherwise each vertex of X n Bq is isolated in G[{ai, . . . , a\x\} U X]. Consequently 
there are at most 2 fc ~Lfc/M)J < 2 1 +( 1 - 1 /M))fc scts X with t[X] > of cardinality at most (1 - l/a)k. At 
the same time there are at most k{rf./ a -i) sets of cardinality greater than (1 — l/a)k. By using the binary 
entropy function, we get (r A / a n) = £>*(2 H(1/Q)fe ), where H{p) = -plog 2 p- (1 -p) log 2 (l— p). For d > 2 and 
a = 3.55 we have 2 ff ( 1 / Q ' < 2 X ~ Consequently if we skip the computation of values i[A] for sets X of 
cardinality at most (1 — l/a)k, such that J H 5o ^ 6, we obtain the claimed running time, which finishes 
the proof of Theorem 11.31 

Note that the constant a = 3.55 can be improved if we have a stronger lower bound on d. However, in 
our analysis it is crucial that a > 2. 

7 Conclusions and open problems 

We would like to conclude with two open problems that arise from our work. First, can our ideas be applied 
to obtain an 0*(2^~ e ^ n ) time algorithm for computing the chromatic number of graphs of bounded average 
degree? For graphs of bounded maximum degree such an algorithm is due to Bjorklund et al. 0]. 

Second, can we make a similar improvements as in our work if only polynomial space is allowed? To the 
best of our knowledge, this question remains open even in graphs of bounded maximum degree. 
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